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ABSTRACT 

A  generalized  Morse  index  theory  is  used  to  study  the  existence 
of  travelling  wave  solutions  of  a  diffusion-reaction  system  of 
equations.  The  reaction  system  is  assumed  to  be  "close"  to  one 
which  admits  an  attracting  manifold  of  critical  points.  A  scaling 
argument  is  used  to  see  that  the  equations  for  travelling  waves  of 
the  full  system  are  then  close  to  a  system  with  a  normally  hyper¬ 
bolic  manifold  of  critical  points. 

Standard  perturbation  theorems  are  already  available  to  study 
the  behavior  of  solutions  of  the  "perturbed"  system  which  lie  near  the 
critical  manifold  in  terms  of  a  (derived)  system  of  "slow"  equations 
on  the  manifold  itself.  Here,  another  such  theorem,  dealing  with 
aspects  of  the  system  which  can  be  described  in  terms  of  isolated 
invariant  sets,  is  proved.  Specifically,  it  states  that  isolated 
invariant  sets  of  the  slow  equations  correspond  to  isolated  invariant 
sets  of  the  full  system,  and  that  the  Morse  index  of  the  latter  set 
is  an  n-fold  suspension  of  that  of  the  former  where  n  is  the 
number  of  unstable  normal  directions. 

These  theorems  are  applied  to  a  standard  continuous  space- 
time  natural  selection-migration  model  for  a  diploid  organism 
when  the  selective  strength  is  weak.  The  selection  is  assumed  to 
be  determined  by  a  single  locus  at  which  the  number  of  available 
alleles  is  arbitrary,  and  the  critical  manifolds  are  found  in  this 
case. 

The  perturbation  theorem  is  applied  to  a  system  with  only  two 
alleles  in  a  situation  where  the  existence  of  a  travelling  wave  for 
the  slow  equation  has  long  been  known.  The  conclusion  is  that  the 
full  system  also  admits  a  corresponding  travelling  wave.  The  index 
theory  is  of  use  because  the  travelling  wave  itself  is  part  of  an 
isolated  invariant  set. 

AMS  (MOS)  Subject  Classification:  34C05,35B99,  92A17 
Key  Words:  Travelling  Waves,  Diffusion  Reaction,  Population 
Genetics,  Morse  Index,  Critical  Manifolds 
Work  Unit  Number  1  (Applied  Analysis) 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-S0- 
C-0041  and  National  Science  Foundation  Grant  Nos.  MCS800-1816 
(Conley)  and  MCS79-04443  (Fife) . 


The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MRC,  and  not  with  the  authors  of 
this  report. 


CRITICAL  MANIFOLDS,  TRAVELLING  WAVES  AND  AN  EXAMPLE 
FROM  POPULATION  GENETICS 
C.  Conlevy  and  P.  Fife 

§  1.  Introduction . 

A  critical  manifold  for  a  system  of  ordinary  differential  equations  means 
a  manifold  of  critical  points.  The  system  of  equations  from  population  genetics 
which  is  treated  here  admits  such  manifolds  when  no  "selection"  terms  are 
present. 

General  perturbation  theorems  are  available  for  the  study  of  the  equations 
with  small  selection  terms  added.  In  §  2  a  simple  version  of  such  a  theorem, 
which  is  particularly  adapted  to  the  motivating  question  of  this  report,  is  proven. 

The  results  of  this  report  are  applied  to  prove  the  existence  of  travelling 
wave  solutions  of  the  system  of  partial  differential  equations  representing  a 
standard  continuous  space-time  natural  selection-migration  model  for  a  diploid 
organism,  when  the  selective  strength  is  weak.  Selection  operates  at  a  single 
locus  with  two  alleles,  so  the  system  has  three  equations  (one  for  each  genotype); 
however  most  of  our  formalism  applies  to  the  multi-allele  case.  Underdominance, 
when  the  two  homozygotes  are  fitter  than  the  heterozygote,  is  assumed,  since 
that  is  the  case  when  a  unique  travelling  wave  front  is  expected . 

The  key  step  in  the  proof  of  the  result  is  a  scaling  which  presents  the 
travelling  wave  equations  as  a  small  perturbation  of  a  system  with  a  critical 
manifold.  The  appropriate  limiting  problem  on  this  manifold  has  been  well 
studied  and  is  known  to  have  a  unique  solution. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29- 
80-C-0041  and  National  Science  Foundation  Grant  Nos.  MCS800- 
1816  (Conley)  and  MCS79-04443  (Fife)  . 


With  the  aid  of  the  results  of  §  2,  the  main  existence  result  is  proved 
in  §4.  The  proof  can  be  played  back  to  a  shooting  argument  or  an  "index" 
argument.  Since  the  latter  has  been  carried  out  in  [1]  the  details  are 
omitted  here .  However,  a  fairly  detailed  outline  of  the  way  in  which  the  index 
arguments  can  be  used  is  given  even  though  the  simplicity  of  the  example  treated 
doesn't  really  warrant  it.  The  point  is  to  set  things  up  for  later  treatment  of 
more  complicated  situations. 
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§  2.  Behavior  of  solutions  near  a  critical  manifold, 
2.1.  Critical  Manifolds . 

A.  Definition . 

Let 


(1)  du/dt  =  U(u) 

be  a  system  of  differential  equations  with  u  in  ]Rn  .  A  critical  manifold, 
to  ,  for  this  system  means  a  manifold  of  critical  points  ( i.  e.  points  at 
which  U  =  0  ) . 

The  critical  manifold  is  non-degerate  ct  u  e  to  if  the  matrix 
9U/9  u  evaluated  at  u  has,  in  addition  to  the  zero  eigenvalues 
corresponding  to  directions  in  to ,  only  non-zero  eigenvalues. 

If  these  additional  eigenvalues  have  only  non- zero  real  parts,  the 
manifold  is  called  normally  hyperbolic  at  u ;  if  only  negative  real  parts, 
attracting  at  u . 

Though  we  will  not  use  the  fact,  it  is  easy  to  see  that  the  existence  of  an 
invariant  manifold  is  inherited  by  the  tangent  equations: 

B.  Theorem :  Suppose  U  is  C*;  then  to  (1)  there 
corresponds  the  system  of  tangent  equations : 

du/dt  -  U(u) 

(2) 

dv/dt  =  (9  U/3  u) v, 


For  this  system  the  set  of  points  to  -  (  (u,  v)  |  u  c  to  and  v  is 
tangent  to  to  at  u  ]  is  a  critical  manifold  with  twice  the  dimension  of  to  . 


If  to  is  non-degenerate,  hyperbolic  or  attracting  at  Uq  then  the 
same  is  true  of  to  at  points  (uQ,  v)  in  to. 


Proof .  Since  U  |  to  =  0  ,  if  u  e  to  and  v  is  tangent  to  to 

then  U(u)  =  0  and  (9U/3u)v=  0.  Therefore,  fn  is  a  critical 

manifold.  The  analogue  of  3U/3u  in  the  definition  is  the  matrix 


where  the  *  refers  to  terms  irrelevant  to  the  present  issue.  The  eigenvalues  of 
this  matrix  are  the  same  as  those  of  9  U/  9  u  and  so  the  second  statement 
of  the  theorem  is  true. 


C.  The  existence  of  critical  manifolds  for  (1)  also  implies 
existence  of  such  for  the  "projectivized"  tangent  equation  but  a  precise  theorem 
requires  a  more  refined  hypothesis.  Since  it  is  not  needed  here,  it  is  omitted. 

The  reason  for  mentioning  the  tangent  and  projectivized  tangent  equations 
is  discussed  in  2.  3D. 


2.  2  A  Parametrized  Family  of  Equations. 

The  perturbation  theorem  concerns  the  smooth,  one  parameter  family 


of  equations 


u=Ue(u,e)=U0(u)  +  eU(u,e)  (*  = 


and  In  particular  the  behavior  of  solutions  near  a  critical  manifold  which  is 
assumed  to  exist  when  e  =  0  . 

Letting  x  denote  coordinates  in  the  manifold  and  y  complementary 
coordinates,  the  equations  can  be  written  as  : 
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(4) 


x=  X£(x)y  +  Xq(x,  y)  +  eX(x,  y,  e) 
y=  Yt  (x)Y  +  Yq(x,  y)  +  8  Y(x,  y,  e) 


where  X  £  and  Y£  are  matrix  valued  functions  of  x ,  and  Xq  and 
Yq  are  of  order  at  least  two  in  y.  These  terms  come  from  UQ  while 
the  last  come  from  e  U  . 

Now  the  critical  manifold  is  non-degenerate  at  x  if  and  only  if  Y£  (x ) 
is  non- singular.  If  this  non- degeneracy  condition  is  satisfied  for  x  in 
some  compact  subset,  say  N ,  of  the  critical  manifold,  then  coordinates 

can  be  chosen  so  that  X£  ( x )  =  0  on  N . 

Namely,  let  x  =  x  -  X£  (x)  Y£~*y.  Then,  computing  modulo 
terms  of  order  two  in  y  or  one  in  e  ,  there  results : 
x  =  x  -  X£  Y"1  y  =  X£  y  -  X£  Y  Y{  y  =  0  .  The  equations  (4 )  then 
assume  the  form  j 


(5) 


x  =  Xq(x,  y)  +  6  X(x,  y,  e) 

y=  Y£  (x  )y  +Yq(x,y)+eY(x,y,e), 


with  X 

q 


and  Y„ 

q 


of  order  at  least  two  in  y . 


Assume  now  that  the  critical  manifold  is  hyperbolic  for  points  x 
in  a  compact  set  N.  Then  the  y-coordinate  can  be  modified  if  necessary 


so  that  for  x  e  N,  Y£(x)  has  the  form: 


(6) 


Vx>  s 


Z,(x) 


W  £(x)J 


where  Z{  (x)  is  negative  definite  and  W{  (x)  is  positive  definite. 
Let  y  =  ( z  ,  w )  where  the  decomposition  corresponds  to  that  of  Y 


and  let 


i 

W  and  Z ,  W  be  the  matching  notation  for  Y  and 

M  4 


t 


t 


L 


I 


-5- 


Y .  Then  the  equations  have  the  form: 


x  =  Xq(x,  y)  +  e  X(x,  y,  e) 

(7)  z  =  Z£  (x)z  +  Zq(x,  y)  +  e  Z(x,  y,  e) 

w  =  W£(x)w  +  Wq(xy)  +  e W(x,  y ,  s ) 

The  dim  (z)  will  be  called  the  number  of  stable  directions;  dim  ( w )  the 
number  of  unstable  ones, 

2.3  Compact  Invariant  sets  near  the  critical  manifold . 

A.  Theorem:  Suppose  the  critical  manifold  in  Is  hyperbolic 
for  all  points  In  some  compact  set  N .  Then  there  is  a  neighborhood  in 
u- space  of  the  interior  of  N  (rel.  in)  which  is  independent  of  e  and 
such  that  the  distance  from  invariant  sets  in  this  neighborhood  to  points 
in  m  is  of  order  e  . 

More  specifically,  suppose  the  critical  manifold  m  =  {  u  |  u  =  ( x ,  0  ) } 
of  (5  )  is  hyperbolic  at  points  in  the  compact  set  N  c:  m .  Using  the 
coordinates  of  (7),  for  k>  0  define  N(k)  =  (u  =  (x,z,w)| 

1  z  I  »  M  *  M- 

Then  there  are  positive  constants  eQ,  kQ  and  kj  such  that: 
if  e  ^  eQ  and  S  is  an  invariant  set  in  N(kj),  then  S  is  in 
N(ek0). 

In  fact,  if  e  s  e0>  ke  [ekQ,  kj]  and  u  e  N{k),  then 
whenever  |z|  =  k,  |z|*<0  and  whenever  J  w  J  =  k,  J  w  |  *  >  0  • 


Proof  .  The  last  sentence  implies  the  previous  one  as  follows.  Assuming 
S  is  in  N(kj),  let  k  be  the  smallest  number  such  that  SC  N(k). 

Then  (if  S  isn't  empty)  there  must  be  a  point  of  S  at  which  JzJ  =  k 
or  |w|  -  k.  Now  if  k  were  in  [ekg,  k^]  then  either  |z|*  <  0 
or  |  w  |  *  >  0  —  whichever  is  relevant  —  and  it  would  follow  that  the 
solution  leaves  N  ( k )  in  one  or  the  other  time  direction.  Therefore  k 
cannot  be  in  [e  kQ,  k^]  .  Since  k  s  k^,  it  must  be  less  than  ekg  . 

Thus  the  last  sentence  implies  the  full  theorem 

Now  let  Sq  and  k'j  be  positive  constants  such  that  the  right 
hand  sides  of  (7)  are  defined  for  e  ^  eQ  and  |  z  |  ,  |w|  s  kj  . 

Let  Cj  be  a  bound  for  [Z|  and  |W  (  for  all  ue  N(k'1) 
and  e  s  eQ  . 

Choose  kQ  so  that  kg*(2  +  c^)  is  a  lower  bound  for  the  quadratic 
forms  -ZSL  (x)  and  (x)  whenever  x  e  N.  Let  cQ  be  a 
bound  for  jZ  j/]yj2  and  |w  l/|y|2  for  all  u  e  N(k',). 

Define  k^  =  min  (k'^ ,  1/cgkg).  Now  suppose  e  ^  eQ  , 

k  e  [ek0,kj]  ,  u  e  N(k)  and,  say,  Jz|  =  k.  Compute: 

<il*|V  =  (z,Z,  z)  +  <z,Zq)  +  eU.Z)  s  ->=0 1(2  +  c1)  |Z|  2  +  |z|  |Zq|  + 

e  |  z  |  |  Z  |  .  Since  u  e  N  ( k )  and  |  z  [  -  k  ,  jw|  s  k  and 
|  y  ]  2  s  2  k2.  Therefore  : 

( ■■  l  z  I  )  s  I  z  I  t  "^o  (2  +  c^)k  +  2k  Cq  +  ecj]  - 

*  |  z  [  [  2 k  (  -kg  +  k Cg  )  +  Cg  ( -kg  k  +  e  )]  . 
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Now  since  k  e  [  ekQ,  k^]  ,  -kQ"^k+  e  <  0.  Also,  because 
kj  s  1/  c0kQ  ,  -kj1  +  kcQ  <  0  .  Furthermore  at  least  one  of  these  terms 
is  negative,  so  (i  |  z  (  2  [  *  <  0  . 

If  |w|  =  k,  a  similar  computation  shows  that  (i|wj2)*  >  0 
and  the  theorem  follows . 

B.  Definition  .  A  non-degenerate  isolating  block  for  the  differential 

equation . 

(8)  v  =  f(v) 

means  a  compact  set  B  described  in  terms  of  a  finite  number  of  (  smooth  ) 
functions  Li»  •  •  • »  Lm  as  B  =  {v|  L^v)  s  0  for  i  =  1 ,  . . . ,  m )  . 

It  is  also  required  that  if  v  e  9  B ,  then  for  one  of  the  defining  Li  , 

L.(v)  =  0  and  Li  ( v )  (  =  V  L  •  f  (v) )  /  0  .  If  (this)  L^v) 
is  positive  at  v ,  v  is  called  an  exit  point;  if  negative,  v  is  an 
entrance  point.  (It  is  possible  to  be  both  an  exit  and  an  entrance  point.  ) 

Blocks  are  more  generally  defined  in  (1]  ,  page  55,  lines  4  and  5  . 

C.  Theorem:  With  reference  to  (5),  let  Xg  (x  )  =  X(x,  0,0). 
Suppose  B  is  a  non-degenerate  block  for  the  equation 

(9)  x  =  Xg(x)  . 

and  that  fn  is  hyperbolic  at  points  u  =  ( x ,  0  )  with  x  e  B  .  Let 
kQ  be  as  in  Theorem  2.3  A  where  the  N  of  that  theorem  is  replaced 
by  B. 
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Then  there  exists  >  0  such  that  if  e  s  e^,  the  set 

B  =  f  u  =  ( x ,  z ,  w )  |  x  e  B  and  |  z  j  ,  j  w  j  <  e  kQ  }  is  a  block  for  ( 5  ) . 

The  exit  set  of  this  block  is  the  set  of  points  u  =  ( x ,  y ,  z  )  with 

x  in  the  exit  set  of  B  or  |w]  =  ekQ  . 

Proof:  For  the  functions  defining  ¥,  take  the  L^(i  =  l,...,m)  defining 

B,  (but  now  considered  as  functions  of  u  =  (x,  z,  w))  together  with 

Lm+1  =  sk0  “  l2l  and  Lm+2  =  ek0  "  M  * 

Now  given  a  boundary  point  of  B  at  which  [  z  j  =  e  kQ  ,  the  function 
Lnri-1  or  Lnrf2  (resP-  )  is  that  required  in  the  definition.  Namely,  by 
the  last  sentence  of  Theorem  2,3  A  the  derivative  on  the  solution  through 
the  boundary  point  of  the  appropriate  one  of  these  is  not  zero . 

Given  a  boundary  point  of  ]3  with  both  |  z  |  and  |  w  |  /  ekQ  , 
the  x- coordinate  of  this  point  is  in  the  boundary  of  B  .  Let  L,.  be 
the  function  such  that  Lj  ( x )  =  0  and  d  ( Xg  )  ( x )  /  0  . 

The  derivative  of  on  solutions  of  ( 5  )  is  given  by 
dLi(Xq(X’  y)  +  6  X(x»  y.  e))  =  edLi(Xs(x))+  dL^X^fx,  y)  +  e(X(x,  y,  e) 
X(x,  0,0)]  .  Now  in  B,  |y  |  2  s  2  e2  k  q  .  Since  Xq (x,  y)  is  of 
order  at  least  two  in  y  and  X(x,y,e)-X(x,0,0)  is  of  order  at 
least  one  in  y,  the  terms  in  square  brackets  are  of  order  at  least  two  in 
s.  Because  dL^(Xg(x))/  0  the  derivative  of  on  solutions  of 
(5)  is  non- zero  in  a  neighborhood  of  u  =  (x,  y,  z)  provided  e  is 
small  enough.  Covering  the  boundary  of  B  with  a  finite  number  of  such 
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neighborh oods  and  taking  to  be  the  least  of  the  e's  ,  the  result  follows  . 

D.  Remarks :  The  equation  (9)  of  Theorem  2.3  C  above 
provides  (except  for  the  change  in  time  scale )  the  lowest  order  approximation 
for  the  behavior  of  the  x-coordinates  of  solutions  of  (3)  near  in.  Of 
course  these  coordinates  change  at  rate  (at  most)  e  since  m  is  a 
manifold  of  rest  points  of  ( 3  )  when  e  =  0  . 

Definition.  The  equation  x  =  X_(x),  considered  as  an  equation 

s 

on  m ,  will  be  called  the  limit  equation  on  In  ( as  e  —  0  ) , 

Theorem  2.  3  C  states  that  any  properties  of  this  limit  equation  which  can 
be  described  in  terms  of  isolating  blocks(or  "local  Liapounov  functions"  like 
the  Lj )  can  be  carried  over  to  ( 3  ) ,  In  view  of  the  remarks  in  2.  IB  . 
and  C.,  the  same  holds  for  properties  describable  in  terms  of  isolating 
blocks  for  the  tangent  and  projectivized  tangent  equations  and  their  iterates. 

There  are  several  such  properties :  for  example,  the  statement  that 
(9  )  admits  a  compact  hyperbolic  invariant  set  in  a  given  conjugacy  class 
can  be  expressed  in  terms  of  blocks  for  the  tangent  and  projectivised  tangent 
equations.  Thus,  if  S  c  m  is  such  an  invariant  set  of  the  limit  equation, 
then  for  small  e  ,  ( 3  )  admits  a  compact  hyperbolic  invariant  set  conjugate 
to  (and  within  kQ  e  of)  S  .  More  specially,  S  might  be  a  hyperbolic 
critical  point  or  periodic  orbit  of  (9  ),  which  then  perturbs  to  an  orbit  of 
the  same  type  for  ( 3  )  -  though  of  course  the  dimensions  of  the  stable  and 
unstable  manifolds  generally  increase . 


As  another  example,  the  usual  hypothesis  on  a  compact  invariant 
manifold  that  allows  a  perturbation  result  preserving  smoothness  of  order 
r  can  also  be  stated  in  terms  of  isolating  blocks  of  an  iterated  tangent 
equation.  Thus  if  In*  C  to  is  a  compact  invariant  Cr  -  manifold  for  (9) 
which  satisfies  this  hypothesis  then  for  small  e*  (and  assuming  enough 
smoothness)  (3  )  admits  an  invariant  Cr-  manifold  homeomorphic  to  and 
near  m' . 

(]t  is  more  challenging  to  try  to  find  a  perturbable  property  not  describable 
in  terms  of  blocks  .  ) 

In  a  different  direction,  it  is  (  also )  well  known  that  with  a  sequence 
of  coordinate  transformations,  (9  )  can  be  replaced  by  better  approximations 
for  the  behavior  of  the  x-coordinates  of  those  solutions  in  an  invariant  set  S 
near  In .  In  fact  In  itself  can  be  replaced  by  critical  manifolds  of  the 
transformed  equations  which  pass  within  0(en)  of  S.  ( In  other  words, 
uniform  asymptotic  expansions  can  be  given  for  solutions  in  S.  ) 

The  use  of  Theorem  2.  3  C  in  the  setting  of  the  Morse  index  is  illustrated 
§3.4  B  .  In  [3]  a  more  general  theorem  about  isolating  neighborhoods  for 
parametrized  families  of  equations  is  proved.  It  can  be  used  in  connection  with 
the  Morse  index  in  the  same  way  . 


2. 4  Travelling  waves . 


A .  The  travelling  wave  equation . 


The  motivating  problem  of  this  report  concerns  the  system  of  travelling 
wave  equations  associated  to  a  partial  differential  equation  of  the  form : 

(10)  3  u/at  =  a2u/af2  +  U0(u)  +  a2U(u,  a). 

The  travelling  wave  equations  are  ordinary  differential  equations  derived  on 
assuming  a  solution  of  (10)  in  the  form  u(t,  £)  =  u(£-act). 

Letting  t  =  £  -  c  t ,  and  1  =  d  /  d  r  and  dropping  the  tilde, 
these  equations  assume  the  form 

(11)  -ocu'  *  u"  +  UQ  ( u)  +  c2  U(u,  a)  . 

We  are  interested  in  the  case  where,  when  o  =  0  ,  the  equations 
for  the  spatially  independent  solutions  of  (10  )  (i.  e.  the  equations  without 
the  diffusion  terms )  admit  a  non-degenerate  critical  manifold.  As  in  2.2, 
this  implies  there  is  a  coordinate  transformation  u  i — >  (x ,  y )  such 
that  these  equations,  namely 

(12)  u  =  UQ(u)  +  c2  U(u,  a). 


transform  to 

(13) 


x  =  Xq(x,  y)  +  o2  X (x ,  y,  c) 
y  =  (x)  y  +  Yq(x,  y)  +  c2  (x,  y,  c) 


where  X  and  Y  have  order  at  least  two  in  y . 

4  4 
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Now,  in  general,  the  coordinate  transformation  u 


*(z  ) 


transforms  ( 12  )  to : 

(14)  z’  =  (d*)"1  UQ(u(z))  +  oZd9'1  U(u,  o) 

where  d  <p  is  the  Jacobean  matrix  of  the  transformation.  Therefore,  if 
z  =  (x,  y)  and  <p  is  the  transformation  u  — >  (x,  y)  above,  then 
(d<p  )"*  UQ(u  (z  ) )  +  aZ  d<p~^  U(u,  a)  is  the  right  hand  side  of  (13). 

Applying  the  transformation  u  =  tp  ( z  )  to  ( 11 )  ^  one  computes  : 
u*  =  d<p[z  )  , 

u"  =  dZ<p(z' )  +  d*(zM)  and,  from  (11), 

«ocdf(z')  =  d2«p(z’)  +  df(z,,)-f-Uo(u(z))+02U(u(z),  a),  or 
z"  =  -ac  z  -  (d*)"1  UQ  (u(z))  -  aZ(d<p)~1  U(u(z),  a)  -  d<p~l  {dZ  <p  (z  )) . 


Now  let  z=  (x,  ay)  (rather  than  (x,y)).  Then  the  terms 
-d y"1  U0  ( u ( z ) )  -  a2  df"1  U(u(z),  a)  take  the  form  ( from  the  right  hand 
side  of  ( 13 )  with  the  noted  modification ) : 


(16) 


Xq(x,  ay)  +  a2  X(x,  ay,  a) 

Yf  (x)  ay  +  Yq(x,  ay)  +  a2Y(x,  ay,  a) 


In  particular,  excepting  Y.  (x)  ay,  these  terms  are  of  order  at  least  two 

2 

in  a  and,  when  y=  0,  the  x- term  reduces  to  a  X(x,  0,0);  this 

is,  in  the  notation  of  Theorem  2.  3  C  ,  X  (x  ) . 

s 

On  introducing  £  and  tj  ( corresponding  to  x  and  y)  by 
(o  £  ,  o  T  )  =  z  ,  the  last  equation  of  (14)  now  leads  to  the  system: 
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L. 


Nate*** 


wrmmt  »I,  -r. 


(17) 


x'  =  a  f* 
y'  =  n 

=  -oc?  -  a-1Xq(x,  ay)  +  aX(x,  ay,  a)  -  (§  ,  tjj 

V  =  -adl  -Y£(x)y  +  a"1  Yq  (x,  ay )  +  aY(x,  ay ,  a)  -  ( £  ,  rj ) 

where  d^p  (d  <j»  (  a~z' ) )  is  decomposed  into  <E>  ( ^  ,  )  and  ^  T  ) 

in  accordance  with  the  decomposition  a-*z'=  ,  t\  ) .  ( <£  and  ^ 

also  depend  on  x  and  y  of  course ;  they  are  quadratic  in  £  and  tj  .  ) 


B.  Theorem ;  Suppose  the  equations  (12)  for  the  spatially 
independent  solutions  of  (10),  that  is  of  du /  dt  =  d  2u/d  £  2  +  U0  (u)  +  a2  U(u,  a), 

admit  a  non- degenerate  critical  manifold ,  In  ,  when  a  =  0  .  Then  the 
same  is  true  of  the  equations  for  travelling  waves  if  they  are  appropriately 
scaled  ( as  they  are  above ) . 


More  specifically,  when  a  =  0  ,  the  equations  ( 17 )  admit  the  critical 
manifold  ft  =  {(x,y,  £,T])|y  =  T)  =  0}.  Coordinates  on  ft  are 
given  by  (x,  £  )  and  the  limit  equations  on  ft  are  given  by: 


(18) 


x'  =  i 

=  -C|  -  Xs(x)  -  *(! ,  0) , 


where  the  limit  equations  on  ft  are : 
(19)  x'  =  X  s  (x) . 
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The  critical  manifold ,  fiT  ,  is  non-degenerate  if  fn  is,  and  is 
hyperbolic  if  tn  is  attracting.  (If  the  linearized  equations  for  in 
have  a  real  positive  eigenvalue ,  to  is  not  hyperbolic ) .  If  In  has 


n  stable  directions ,  in  has  n  unstable  ones . 


Proof;  On  setting  a  =  y  =  ^  =  0  ,  the  right  hand  side  of  (IT)  vanishes 

( recall  that  X  and  Y  are  of  order  2  in  their  second  variable 
q  q 

so  a-1X  (x,  ay)  and  c-1Y  (x,  ay)  are  of  order  one  in  a).  The 

M  tL 

relevant  portion  of  the  matrix  for  the  linearized  equations  is 


It  is  non-degenerate  if  Yf  (x)  is  and  it's  eigenvalues  are  the  square 
roots  of  those  of  Y(  (x). 

Setting  y  =  =  0  and  taking  the  terms  of  order  a  in  the 

x,  £  equations  give*  (18).  The  equations  (19)  come  on  taking  the 
lowest  order  terms  in  the  corresponding  equations  for  (12  ).  These  are 
of  order  a2  . 


§  3.  Some  Equations  from  Mathematical  Population  Genetics. 

2.1  The  Growth  Equations . 

As  we  consider  continuous -time  models  of  natural  selection  at  a  single 
gene  locus.  Equation  (21)  and  equations  analogous,  in  varying  degrees,  to 
those  in  Sections  3.  IB,  C  can  be  found  in  a  number  of  works;  see,  for  example, 
[4-8]. 

There  is  supposed  to  be  a  population  with  a  total  of  T  individuals  which 
is  classified  according  to  which  two  of  n  possible  "alleles"  appear  at  a  particular 
gene  locus.  The  number  of  alleles  of  type  is  2N^ ,  and  u^(  =  u^  )  is  the 

number  of  ordered  A^A^  genotypes.  Consequently,  and 

T  5  2iNi  *  suuir 

The  growth  equations  are  assumed  to  have  the  form: 

(2!)  «ij  =  /ST^NjNj-  6  1 AAj  +  eg^. 

where  p  and  6  can  be  taken  to  be  ( positive )  functions  of  u  and 
g  i  j  is  some  function  of  u  ( frequently  a  constant  multiple  of  u  i  ^ ) . 

In  particular,  the  birth  rate ,  0T"1  ,  of  the  u^  types  is  proportional 

to  the  product  of  the  numbers  of  individuals  with  the  i^  and  1th  alleles 
(respectively).  This  is  the  "random  mating"  hypothesis.  The  (essential) 
restriction  that  p  and  the  exponential  death  rate  ,  6  ,  do  not  vary 

from  equation  to  equation  corresponds  to  an  assumption  that  differences  in 
fertility  and  survival  ability,  accounted  for  in  the  gA j  terms,  are  small 
enough  to  be  treated  as  a  perturbation,  hence,  the  parameter  e  .  These  same 
terms  may  account  for  small  deviations  from  random  mating  and  for  the  systematic 
introduction  into  the  population,  through  mutation  or  otherwise,  of  small  numbers 
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of  the  various  genotypes  . 


B .  Random  Mating  and  the  Hardy-Weinberg  Surface  . 

In  view  of  the  random  mating  hypothesis  it  is  easy  to  guess  that 
(if  e  =  0)  the  population  might  tend  to  a  state  where  =  T~1N1 

(the  alleles  should  be  independently  distributed) .  This  suggests  introducing 
the  functions 

(22)  hij  =  T"1  Ni  Nj  -  u^  . 

Writing  u  for  u^  etc.  and  p  for  0-6  (the  effective 
exponential  growth  rate )  the  equations  become : 

(23)  u=pu+0h  +  eg. 

The  derivative  of  h  along  solutions  of  this  equation  is 

• 

h  =  pdh(u)  +  0dh(h)  +  edh(g).  Since  h  is  homogeneous  of 

degree  1,  the  first  term  on  the  right  is  just  ph  .  Also,  since  Zjh^  =  0, 

Nj  =  2^  uij  and  T  =  Z^  do  not  change  under  displacements  in  the 

h  direction;  therefore  dh(h)  =  dh(-u1j)  =  -h. 

Recalling  that  p  -  0  =  -  6  ,  this  gives  : 

(24)  h  =  -6h  +  eah(g). 

Of  course  the  death  rate  ,  5  ,  is  positive,  so  if  e  -  0  h  tends  to 
zero  exponentially . 

The  surface  defined  by  h  =  0  is  called  here  the  Hardy-Weinberg  surface. 
It  is  n-dimensional  and  coordinates  on  it  are  supplied  by  the  vector 
N=  (Nj,...,Nn).  In  fact,  if  fi  denotes  the  set  of  ( n2  -  n )  h  ^  j 


with  i  /  j ,  then  when  T  /  0  the  transformation  u 
well  defined  and  invertible  and  has  Jacobean  one . 


(N,  h  )  is 


Namely,  from  the  definition  (22),  liij  =  T-1  Nj.  Nj  -  uIj  '  ** 


follows  that  Ujj=T  N^Nj 


so  that  N  and  h  determine 


the  with  i  /  j  .  Then  from  one  can  determine 


To  compute  the  Jacobean,  observe  that  dNj  A...  A  dNfl  Adh^ 

must  be  the  same  as  dN,  A...  A  dN  Adu..  since  h..  and 

i  n  ij  ij 

u1;j  differ  by  a  function  of  the  N' s  .  Then  from  Nj  =  2j  uij  At 
follows  that  dNjA...  A  dNfi  A  (/V^jdu^)  A,  ^  d  u^  j  since 
the  terms  in  Nj  of  the  form  d  u  j  j  with  i  /  j  already  appear  in 
the  product  AiA  dUjj  .  The  results  follows. 


C  .  Critical  Manifolds . 

On  this  Hardy- Weinberg  surface ,  the  equations  ( with  e  =  0 )  are 
just  u  =  p  u  .  Therefore  the  zero'  s  of  the  effective  growth  rate,  p  , 
provide  critical  points.  The  growth  rate  is  obviously  positive  for  some 
values  of  u ,  and  is  negative  if  u  is  so  large  that  the  environment 
can  not  support  the  total  population.  Therefore  the  Hardy-Weinberg  surface  must 
be  separated  by  a  set  of  points  where  p  =  0  .  This  set  is  assumed  to  be  a 
manifold,  and  is  therefore  an  n-1  dinensional  one . 
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It  has  already  been  seen  ( from  ( 24 ) )  that  ( when  e  =  0  ) ,  this 
manifold  is  attracting  in  the  n  (n  -  l)/2  independent  directions  corresponding 
to  the  h  ^  with  i  ^  j  .  The  situation  in  the  remaining  direction  depends 
on  the  derivative  of  p  . 

An  easy  case  to  treat  is  that  where  p  depends  only  on  the  total 
population ,  T ,  and  has  a  zero  at  some  critical  value  TQ  .  Then  part 
of  the  critical  manifold,  Hi ,  is  the  set  {ujhjj  =  0;  T=T0}. 

Coordinates  in  this  manifold  are  supplied  by  the  gene  frequencies : 

(25)  P*  =  Nj/T  . 

Now  define  coordinates  P,  and  t  by: 

P  =  (Px,  Pn) 

(26) 

t  =  T  “  T0 

and  define  ,  G ,  the  constant  X  and  the  function  r  ( t ) 

( of  order  two  in  t  as  t  —  0  )  by : 


(27) 

*  =  -[p(T0  +  t)(T0  +  t)]  (T0)  and 
r(t)  =  p(T0  +  t)(T0  +  t)  +  Xt  . 


Then  the  equations  ( 21 )  assume  the  form : 


?i  =  e(TQ  +  t)-i[  Gj  -  G  P.] 

(28  )  t  =  -X  t  +  r (t)  +  e  G 

h  =  -6h  +  edh(g). 

To  see  this  add  (21)  over  j  to  get  =  pN+eG^  and  then  over  i 

to  get  T  -  pT  +  e  G  .  Then  P  =  (T-1  N)‘=  eT'1^  -  GPj]  as 

required.  The  equation  for  t  comes  from  (26)  and  (27)  with 

•  • 

t  -  T=p(T)T  +  eG.  That  for  h  comes  from  (24).  The  equations 
for  h  are  redundant  of  course  and  h  should  be  interpreted  as  the 
hAj  with  i  /  j,  i.  e.  “’ll  .  Observe  that  the  equations  (28)  have 
the  form  of  ( 5 )  with  P  replacing  x  and  ( t ,  £ )  replacing  y  . 

The  analogue  of  Yf  is  a  diagonal  matrix  with  a  -  X  on  the  diagonal 
and  remaining  diagonal  entries  equal  to  -6  .  The  analogue  of  X  is 

'“I 

zero  and  that  of  Yq  is  r(t)  . 

It  will  be  assumed  that  \  >  0  so  this  critical  manifold  is  attracting. 

The  limit  equations  on  the  critical  manifold  tn=[P,t,1i[t=fr=0} 
are  given  by : 

(29)  ^  =  (G.  -  GP.)/T0  . 

With  appropriate  choice  of  g^  ,  these  are  the  well-known  Fisher  equations  . 

D .  The  Travelling  wave  equation . 

The  travelling  wave  equations  associated  to  (23  )  are  (defining  v 
by  the  equation  ) : 


-20- 


U  =  V 

(30) 

v'  =  -0cv  -  pu  -  0  h  -  a2  g , 

2 

where  e  has  been  replaced  by  a  and  the  wave  velocity  is  equal  to 
ac . 

The  equations  in  the  coordinates  (P,  t,^)  and  their  derivatives 
could  be  derived  from  the  general  form  (17)  or,  since  they  are  the  only  ones 
required,  those  for  P  and  Q=  a"*  p'  can  be  computed  directly  as  follows. 


Defining 

t'  =  T  =  Z^v^  to  be 

a  w 

there 

follows 

«  i  2 

aw  =  vij  =  -  o’  cw 

-  pT  - 

a2  G  . 

'  i 

Nt  =  (TPi)  =  awPi  +  aTQj.  It  then  follows  that 

Nj'  =  aw'  Pj  +  2oZwQ.  +  aT .  Rearranging  and  using  the  above 
expression  for  aw  gives  QA  =  (  aT )"  [N'^  +  a  cwPi  +  pTPi  + 
a^Pj  -  2a2wQ.  ]  . 

Now  using  Sj  v^  =  n'^  and  adding  the  second  equation  of  (30  )  over 
j  gives  Nj"  =  >acNj  -  pTPj  -  a2Gi.  Using  N^  =  awpj  +  aTQ1 
there  results  N^'  =  ■a2cwP1  +  a2cTQi  -  pTPj  -  a2Gi. 

Substituting  this  last  expression  for  N1^  into  the  preceding  equation 
for  ,  there  follows :  =  (  aT  )"*  [ -  a2  cT -  a2  Gi  +  a2GPj  -  2  oLwQ.]  . 

Thus  the  equations  for  P  and  Q  are : 

=  oQt 

(31)  .  _i  ,  _i 

Qi  =  ~ocQi  -  o T  X[G.  -  GPj]  -  ZaT  1wQi. 
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These  are  the  analogues  of  the  equations  for  x  and  £  in  (18). 
However,  it  should  be  noted  that  there  should  really  be  only  2  n  -  2 
equations  in  (31)  since  siPi  ~  1  and  consequently  ~  0  • 

The  critical  manifold  for  the  travelling  wave  equations  (23  )  with 
e  -  0  is  the  manifold  of  points  where  t  =  0  (T  =  Tg  ),  w  =  0  and 
h  =  h'  =  0  .  The  limiting  equations  on  this  manifold  are  obtained  by  setting 
these  variables  equal  to  zero  in  ( 31 )  and  taking  the  first  order  terms  in  a  . 

This  gives : 

P'i  =  Q4 

(32) 

Qx  =  -cQj  -  (GA  -  GPp/Tg 

Now  the  critical  manifold  for  (28)  is  obviously  attracting  since  X  >  0 
with  n(n-l)/2  +  1  stable  directions  assuming  u„  =  u^.  By  theorem 
2.4  B  the  corresponding  one  for  (30)  is  hyperbolic  with  n(n  -  1  )/2  +  1 
unstable  directions.  Therefore  (by  2.  3  C)  blocks  for  (32)  carryover 
to  blocks  for  (30). 

In  the  following  example,  this  will  be  used  to  prove  the  existence  of  a 
travelling  wave  for  (30)  from  the  (known)  existence  of  such  a  wave  for  (32), 

3.  2  An  Example  . 

1  2 

A.  The  easiest  example  is  that  when  n  =  2  and  g„  =  u^  +  6^  T  P.P^ 
(the  first  term  representing  a  perturbation  of  the  death  rate,  and  the  second  a 
perturbation  to  the  birth  rate).  Since  Pj  =  1  -  P2  and  Qj  =  -Q2,  only 
one  pair  of  the  equations  in  (3.2)  is  needed. 
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On  the  Hardy -Weinberg  surface  ,  u^  =  TP^  and  on  the  critical  manifold 

1  2 

T  =  TQ  .  Therefore  setting  6^  =  6^  +  6^  and  taking  6^2  =  S21  '  we  ^ave 

(Gj-GP^/Tq  =  P1[61  P1  612  P2  -  611P1  -  2  612P1  P2  -  622p2]  * 

This  last  expression  is  a  function  of  P^  alone,  since  P2  =  1  -  Pj  .  It 
vanishes  when  Pj  =  0  or  1,  and  when  the  6^  are  constant,  it  is  cubic 
with  third  root 

3  =  ^22  '  612  ^  611  +  2612+  622^  * 

The  quantities  1  +  5^  are  called  the  fitnesses  of  the  three  genotypes.  The 
case  we  shall  be  especially  interested  in  is  the  "underdominant"  case 

611  >  612'  622  >  612  * 

Normalizing  the  5^  so  that  §l2  =  0  and  dropping  the  subscript,  (31) 
(for  P^  Q^)  takes  the  form 


(33) 


d  P/ d  t  -  Q 

dQ/dT  =  -cQ  -  P(1  -  P)(P  -  a) 


with  A  =  622/(6u  +  622)  . 

B.  These  equations  have  been  treated  by  several  authors  (e.  g.  [9-13]  ). 
Our  aim  is  to  use  Theorem  2.  3C  to  "lift"  the  results  to  the  corresponding 
equation  for  the  u^  ,  or,  more  specifically,  to  prove : 


Theorem.  Let  u  =  (  uy  )  and  v  =  (vy  )  with  i,  j  =  1,  2  . 

Let  Nj  =  Zj  u  and  T  =  sij  u  ij  •  p  =  p(T)  and 

0  -  0  (T )  be  differentiable  functions  of  T  with  0  and  p  -  p 
both  positive  when  T  is  positive.  Let  h  =  (h^)  =  (T~*  -  u^  ) 


] 


and  let  g  =  g^  be  functions  of  u  which  reduce  to  the  form  6^(T  )  u.^ 
when  h  =  0.  Let  a  and  p  de  real  parameters. 

Assume  that  for  some  TQ  >  0  ,  p(TQ)  =  0  and  [dp/dT]  (TQ)  <  0  . 

Also  assume  that  612(TQ)  =  8  21^0^  =  0  anc*  611^0^  >  822^T0^  >  0 
(under  dominance).  Then  there  is  a  >  0  such  that  to  each  cr  <  ctq 

there  corresponds  a  c  <  0  such  that  if  c  =  c  (a),  the  equations 


du/d  x  =  v 

(34)  2 

dv/dx  =  - ctcv  -  pu  -  ph  -  a  g 

admit  a  non-constant  solution  ( u  p(  x ) ,  v^x))  which  tends  to 

distinct  rest  points  as  x  tends  to  minus  or  plus  infinity. 

On  this  solution  T  =  TQ  +  O  (  a  )  and  the  limits  of  u  at  minus 
and  plus  infinity  ( respectively )  are  given  to  order  a  by : 


(35) 


0  T, 


With  P  =  T~A  (u^  +  Uj2 )  the  function  P^fx)  corresponding 
to  the  above  solution  converges  as  a  —  0  to  a  solution  of  the  equation 
(33)  which  tends  to  the  rest  points  (0,0)  and  (1,0)  as  x  tends 

to  minus  or  plus  infinity  respectively.  (Of  course  c*  (  cr)  must  also 
converge  to  a  value  of  c  for  which  such  a  solution  exists ) . 


Proof:  The  proof  of  this  theorem  could  be  phrased  in  terms  of  the  "shooting' 

method"  since  it  is  quite  easy  to  construct  the  necessary  blocks  explicitly. 


1 


In  fact,  referring  to  ( 33  )  and  the  function  H  given  in  3.  3  A, 
for  all  negative  values  of  c  the  set 

B  =  {  (P,  Q )  1  -  e  ssPsl+e;  H(a,0)+e*H(P,  Q)*H(l,0)+e} 
is  a  block;  it  corresponds  to  the  shaded  region  in  Figure  1.  It  is  easily 

seen  that  for  c  close  to  zero  the 
branch  of  the  unstable  manifold  of  the  rest 
point  (0,0)  leaves  the  block  as 
indicated  by  the  solid  curve,  and  for 
c  large  it  leaves  as  indicated  by  the 
dashed  curve.  It  follows  that  for  some 
value  of  c  it  goes  to  the  rest  point  (1,  0  ) . 

this  block  corresponds  to  one  for  (3.4) 
might  be  constructed. 

Instead,  the  proof  here  will  be  based  on  the  general  facts  about  the 
Morse  index  of  isolated  invariant  sets ;  this  is  to  illustrate  the  use  of  this  machinery 
and  to  set  things  up  for  a  treatment  of  more  complicated  examples.  Of  course  the  same 
features  enter  the  proof,  but  the  index  approach  is  probably  easier  to  implement  in 
general,  particularly  when  it  is  difficult  to  construct  blocks  (or  Wazewski  sets, 
or  local  Liapounov  functions.. , )  explicitly.  On  the  other  hand,  the  index  argument 
rests  on  a  fairly  large  body  of  abstract  definitions  and  theorems .  The  following 
description  of  which  (without  proofs)  is  already  longer  than  a  direct  proof  for 
this  example  would  be .  In  3.  3  A,  B  and  C  and  3. 4  B  and  C  remarks 
pertinent  to  the  example  are  more.  In  3.5  a  general  theorem  for  travelling 
wave  equations  is  given.  Finally  in  3.6  it  is  shown  how  the  machinery  works 
in  the  example. 


0 


Fig.  1 

Now  using  Theorem  2.3  C, 
and  a  similar  "shooting  argument" 
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3.3  Facts  about  (33). 

A.  Lemma  If  c  <  0  the  function  H(P,  Q)=  Q2/2  - 
4  3  2 

P/4  +  (l+a)P/3  -  aP  /  2  is  strictly  increasing  on  non-constant 
solutions  of  (33  ).  Consequently  every  bounded  solution  tends  to  a  rest 
point  in  forward  and  backward  time  and  the  two  rest  points  so  determined 
are  distinct. 

There  are  three  rest  points  for  the  equation;  namely,  ( a ,  0  ) ,  (0,0) 
and  (1,0).  The  values  of  H  at  these  three  points  are  (respectively) 

3 

a  (a  -  2  )/ 12  ,  0  and  (1  -  2a  )/l2  and  increase  strictly  in  the  order 

given  when  0  <  a  <  1/2  (as  is  assumed ) . 

The  set  of  bounded  solutions  is  compact  and  in  fact  for  all  c  is 
contained  in  the  set  ((P,Q)|0  s  P  s  1;  a3(a  -  2  )/12  <  H  <  (1  -  2a  )/l2  }  . 

Proof;  In  terms  of  H ,  the  second  equation  of  (33)  is 
dQ/d  t  =  -cQ  -  9H/  a  P .  This  makes  it  easy  to  see  that,  on  solution, 
dH/dr  =  -cQ  and  from  there  that  H  is  strictly  decreasing  on  non¬ 
constant  solutions.  The  statement  concerning  the  asymptotic  behavior  of  the 
solutions  now  follows  from  well-known  general  facts  about  gradient- like  equations 
(e.g.  see  [1]  ,  Ch.  I,  §6  p.  12  etseq.). 
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Since  the  equation  is  gradient-like  with  respect  to  H,  the  easily 
computed  values  of  H  at  the  critical  points  provide  bounds  for  H  on  the 
set  of  bounded  solutions.  The  bounds  on  P  come  from  the  fact  that  critical 
points  of  P  {  t  )  are  strict  maxima  if  p  <  0  and  strict  minima  if  p  >  1 
(cf.  [1],  Ch.  n  §4. 3  C  p.  31  for  a  more  general  result ) . 

B .  Lemma  There  is  a  constant  k  such  that  if  |  c  J  as  k 
then  the  rest  point  (1,0 )  is  a  component  of  the  set  of  bounded  solutions 
of  (35). 

Proof:  The  point  is  obviously  a  component  of  the  set  of  bounded  solutions 

when  c  =  0  (when  H  is  constant  on  solutions  and  the  set  of  bounded 
solutions  is  easily  pictured ) .  Because  it  is  compact,  the  set  of  bounded 
solutions  depends  upper  semi-continuously  on  c.  Because  (1,0)  is 
a  hyperbolic  rest  point  no  new  bounded  solutions  are  created  near  this  point 
when  c  changes  a  little.  Therefore,  it  remains  a  component  for  JcJ  small. 

C .  Lemma  The  equation 

(36  )  P1  =  P (1  -  P)(P  -  a) 

is  gradient-like  with  respect  to  the  function  P4/4  -  (1  +  a )  P3/3  +  aP2/  2 
and  has  critical  points  at  0 ,  a  and  1  -  a  . 

For  0  <  e  <  a,  the  interval  [a-e,l+e]  is  an  isolating 
block  with  exit  set  (  a  -  e  }  .  It  contains  the  two  rest  points  a  and  1 . 


Proof:  Trivial 


3. 4  Interpretation  of  the  previous  Lemmas  as  statements  about  Morse 
decompositions . 

A .  To  conclude  the  proof  the  notions  of  attractors,  Morse 
decompositions,  the  homotopy  index  and  connection  maps  are  used .  An  attractor 
relative  to  a  compact  invariant  set ,  S ,  means  a  compact  invariant  subset 
which  is  the  go  -limit  set  of  a  neighborhood  of  itself  (Definition  5.1,  p.  32 

Ch.  n,  [1]).  A  Morse  decomposition,  D  =  (  Mj,  M2  ,  . . . ,  Mn  )  ,  of 
S  is  derived  from  an  increasing  sequence  ^  =  Aj,  ...,An  =  S  of 
attractors  relative  to  S  ;  namely,  Mj  is  the  maximal  compact  invariant 
set  in  AjXAj  ^  (Definition  7.1,  p.  40,  Ch.  n,  [1]).  To  each  Morse 
set  Mj  there  corresponds  an  index,  h  ^  which  takes  the  form  of  the 
homotopy  type  of  a  pointed  space  (Definition  5.1  and  5.2,  pp.  51-52, 

Ch.  HI,  [1]  ).  To  a  pair  of  adjacent  Morse  sets,  M  . ,  M  in  a 
decomposition  there  corresponds  a  connection  map,  c  from  h  ^  to  the 
suspension  of  h  ^  ( 7.  2,  7.3  pp.  61-62,  Ch.  m,  [1]  ).  All  of  these 
objects  are  "stable"  under  small  perturbation  (5.3  C,  p.  35,  Ch.  II, 

Thm.  1.4  p.  67,  Thm.  2.5,  p.  70  and  Thm.  3.1,  p.  72,  Ch.  IV,  [1]  ). 

The  application  to  the  present  problem  goes  as  follows  . 

B .  Lemma  By  Lemma  3.  3  A  with  7. 1  C ,  p.  40,  Ch.  II  of 

Cl].  D=  {(1,0),  (0,0),  (a,0)}  is  a  Morse  decomposition  of  the  set 

of  bounded  solutions  of  (  33  )  for  all  values  of  c . 

Furthermore,  if  there  is  no  solution  running  from  (0,  0)  to  (1,0) 

D  =  {(0,  0),  (1,0),  (a, 0)}  is  also  a  Morse  decomposition 


(  7. 1  C  of  [13  again ).  In  particular,  Lemma  3.3  B  says  no  bounded 
solution  runs  to  (1,0),  so  D  is  a  Morse  decomposition  for  small 
negative  c  . 

Also,  Lemma  3.3B  with  Theorem  7.2  C,  p.  &2  of  [1]  implies  the 
connection  map  from  h(  ( a ,  0 ))  to  the  suspension  of  h((l,0))  is  trivial . 

C.  Lemma  By  Lemma  3.3  C,  D1  ={0,1,  a}  is  a  Morse 
decomposition  of  the  set  of  bounded  solutions  of  ( 36  ) .  Since 
[a-e,  1+e]  can  be  deformed  to  its  exit  set  {  a  -  e  )  it  has  index 

o  (cf.  [1],  bottom  of  p.  6  and  top  of  p.  7).  From  (1),  7.  2  B 

p.  62  Ch.  IE  of  [1]  ,  it  then  follows  that  the  connection  map  from  h  (a  ) 
to  h(l)  is  an  isomorphism.  Since  h(l)  =  1  (cf.  middle  p.  6, 

Ch.  I  of  [1]  )  this  connection  map  cannot  be  trivial. 

3. 5  Three  more  abstract  theorems. 

The  idea  now  is  to  "lift"  the  Morse  decompositions  of  3.4  A  and 
3.4  B  up  to  the  equation  (34)  using  Theorem  2.3  C  or  rather  its  following 
consequence  (see  Definition  3.2,  Ch.  HI  p.  45  of  [1]  for  the  definition 
of  isolated  invariant  set  and  isolating  neighborhood,  and  5. 2 ,  Ch.  Ill,  p.  51 
for  the  index ) . 

A.  Theorem ;  Let  to  be  a  critical  manifold  for  the  equation 
(37)  u  =  UQ(u,  e)  +  e  U(u,  e)  (*= 
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when  e  =  0 ,  and  let  S  be  an  Isolated  invariant  set  of  the  limit 
equation  on  in  such  that  in  is  hyperbolic  at  points  of  S .  Then  any 
neighborhood  of  S  in  u- space  contains  a  compact  neighborhood  N 
such  that  for  small  enough  e  ,  N  is  an  isolating  neighborhood  for  (  37  ). 
The  corresponding  family  of  isolated  invariant  sets  will  be  denoted  S  (e). 

If  S  is  contractible  in  a  subset  of  tn  consisting  of  points  at 
which  in  is  hyperbolic,  then  the  index  of  sets  in  N  is  the  product  of 

that  of  S  with  a  pointed  r- sphere  where  r  is  the  number  of  eigenvalues 

with  positive  real  part  of  the  matrix  9  U/  3  u  evaluated  at  any  point  of  S , 
or  in  other  words  the  number  of  unstable  directions  of  Hi  in  N . 

The  family  S  ( e ) ,  will  be  called  an  r-fold  extension  of  S . 

Proof;  This  could  be  played  back  to  Theorem  2.  3C. ,  or  argued  directly  as 
follows.  Let  U  be  a  neighborhood  of  S  in  u-space.  Since  in  is 
hyperbolic  at  points  of  S ,  there  is  an  isolating  neighborhood  N ,  of  S 

( as  an  isolated  invariant  set  of  the  limit  equation )  in  m  H  L  such  that 
m  is  hyperbolic  at  points  in  N .  Then  for  small  enough  k ,  the  set 
N(k)  of  Theorem  2.3  A  is  contained  in  U.  For  small  enough  e  it 
can  be  shown  to  be  an  isolating  neighborhood  for  (37)  as  follows.  Suppose 
given  a  solution  which  stays  in  N (k )  for  all  time  ;  the  thing  to  be  shown 
is  that  this  solution  does  not  pass  through  a  boundary  point  of  N (k ) .  Now 
by  Theorem  2. 3  A,  it  lies  within  0(e)  of  In.  Therefore,  if  e  is 
small,  it  closely  follows  solutions  of  the  limit  equation.  Also,  if  e  is 


small  enough,  the  y  coordinate  of  the  solution  is  smaller  than  k . 

Therefore  if  the  solution  were  to  pass  through  a  boundary  point  of  N(k ) 
it  would  have  to  be  one  whose  x-coordinate  is  in  the  boundary  of  N.  Then 
the  corresponding  solution  of  the  limit  equation  leaves  N  (because  N  is 
an  isolating  neighborhood  of  the  limit  equation ),  But  then  the  given  solution 
would  have  to  follow  it  out  of  N (k )  which  is  a  contradiction.  This 
proves  N(k)  is  an  isolating  neighborhood  of  (37)  if  e  is  small 
enough  (in  [3]  there  is  a  substantial  generalization  of  this  argument). 

To  compute  the  index  of  S  ( e ) ,  a  block  for  S  is  chosen  and  then, 
by  Theorem  2.3  C  a  block  for  S(e)  (e  small).  If  S  is  contractible 

in  a  subset  of  in  at  points  in  which  to  is  hyperbolic,  then  the  block 
for  S  ( e  )  can  be  deformed  to  the  product  of  the  given  one  for  S  and  a 
block  for  a  hyperbolic  rest  point  whose  unstable  manifold  is  the  number  of 
eigenvalues  with  positive  real  part  of  d  U/3  u  evaluated  at  a  point  of  S. 
This  implies  the  statement  about  the  index  of  Sg  (cf.  §4,  Ch.  I  of  [1]  ). 

B.  Theorem :  With  the  hypothesis  of  Theorem  3. 5  A,  suppose 
D  -  {  Mj ,  . . . ,  Mn  )  is  a  Morse  decomposition  of  S .  Then 
D(e)=  {  Mj ( e ) ,  . . . ,  Mn ( s )  )  with  the  ( e )  the  r-fold  extensions 
of  Mj,  is  a  Morse  decomposition  of  S  (e)  for  small  enough  e  . 
Furthermore  the  connection  maps  from  h  ( M  ( e )  to  the  suspension  of 
h(Mj(e))  are  the  r-fold  suspensions  of  the  corresponding  maps  for  D. 
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Proof:  This  is  a  straightforward  consequence  of  the  definitions  and  the 

preceding  theorem. 

The  next  theorem  is  the  result  initiating  §3,  p.  81  of  [2]  It  is 
included  here  for  completeness  and  because  it  can  be  played  back  to  critical 
manifolds  as  follows . 

C  .  Theorem  Let  D  =  {  Mx ,  . . . ,  }  be  a  Morse  decomposition 

of  an  isolated  invariant  set  S  of  the  equation 

(38)  u  =  f (u)  . 

where  u  e  Fr.  Then  for  large  c ,  S  and  D  admit  r-fold  extensions 
to  the  equation 

t 

U  -  V 

(39) 

V  =  -cv  -  f ( u ) 

which  are  obtained  by  realizing  (38)  as  the  limit  equation  for  a  critical 
manifold  of  a  transformed  version  of  ( 39  ),  namely : 

u'  =  e  f (u )  +  ew 

(40) 

w  =  w  -  ef  (u)w  -  e  f  (u)f(u). 

This  critical  manifold  has  dim  ( u )  unstable  directions  . 

Proof:  With  e  =  0,  the  equations  (40)  become  u'  -  0,  w'  =-w 

and  so  admit  Ih  =  {(u,  w)|  w  =  0}  as  a  critical  manifold.  The  limit 

'  The  function  f  ( v )  introduced  on  page  80  of  [  2  ]  can  be  taken  to  be 
v  ( 1  -  v )  ( v  -  a )  with  0  <  a  <  }  . 


equation  is  u'  =  f(u).  Therefore  3.5  A  and  B  apply. 

-2 

The  equation  ( 40 )  transforms  to  39  on  setting  e  =  c  ,  changing 
the  time  scale  by  multiplying  the  right-hand  sides  by  -c,  and  setting 
w  =  -cv  -  f  (u)  . 

3.6  The  conclusion  of  the  proof . 

Using  Theorem  3.5  B  and  Lemma  3.4  B  the  Morse  decomposition 
D  =  {(1,0),  (0,0),  (a, 0))  of  the  set  S  of  bounded  solutions  of 
(33  )  extends  to  a  decomposition  15  (  a)  of  the  extension  S  (  a)  of  S 
to  (34).  It  is  an  elementary  fact  that  the  extensions  of  the  (non-degenerate) 
rest  points  (1 ,  0 )  (0  ,  0  )  and  (a ,  0 )  are  again  rest  points  ( as  well  as 
Morse  sets) .  These  Morse  sets  will  be  denoted  respectively  by  M^(  a), 

M_  (  a)  and  M_  (  a)  . 

U  a 

By  7. 1  C,  p.  40,  Ch.  n  of  [1]  ,  so  long  as  there  is  no  solution 
running  from  MQ  ( a)  to  Mj(a),  the  ordering  of  MQ(a)  and 
Mj(a),  can  be  interchanged  so  that  D=  {  MQ  (  a),  Mj(a),  Ma(a)}  is 
also  a  Morse  decomposition  of  S(  a).  Theorem  3.5  B  with  Lemma  3.4  B 
then  implies  that,  if  c  is  small,  the  connection  map  for  the  pair  Mj  (  a)  , 
M=  ( cr)  is  trivial . 

a 

Now  Lemma  3.  4  C  with  Theorem  3.5  C  and  then  3.5  B  implies  that 
for  large  c ,  D  is  again  a  Morse  decomposition  but  the  connection  map 
for  the  pair  M,  (  a),  M  (  a)  is  no  longer  trivial :  it  is  an  isomorphism 
between  the  index  of  M  (  a)  and  the  suspension  of  the  index  of  M,  (  o) . 

3  A 
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The  latter  is  the  product  of  1  with  a  pointed  one  sphere  (Theorem  3. 5  C ) 
and  then  a  pointed  two  sphere  (Theorem  3.5  B;  cf.  Theorem  2.4  B)  the 
result  is  a  pointed  three  sphere,  so  the  isomorphism  isn't  trivial. 

Now  Theorem  3.1,  p.  72,  Ch.  IV  of  [1]  implies  that 
D  (  a)  =  f  MQ  (  a),  (  a) ,  Ma  (  a) }  cannot  be  a  decomposition  for  all 

c;  in  view  of  7. 1C,  p.  40,  Ch.  HI  of  [1],  this  means  there  is  some 
value  of  c  for  which  there  is  a  solution  connecting  MQ  (  a)  and  M^(  a), 
and  this  completes  the  proof  of  Theorem  3. 2  B  except  for  the  limit  statements 
which  are  obvious . 

3.7  Concluding  remarks . 

The  equations  from  population  genetics  seem  to  be  rich  in  attracting  or 
hyperbolic  critical  manifolds:  in  the  case  of  two  loci  the  "linkage  disequilibrium 
surface"  is  such  a  manifold  if  there  is  a  positive  probability  of  crossovers 
(and  no  selection).  With  several  loci  analogous  statements  are  true  ([4]  ). 

In  the  simplest  case  of  two  loci  and  two  alleles,  the  critical  manifold  of 
the  "spatially  independent"  equations  is  two  dimensional  while  that  for  the 
travelling  waves  is  four  dimensional.  Therefore  the  analysis  of  an  example 
is  more  difficult. 

It  is  not  Just  the  increase  of  dimension  that  creates  difficulty,  though. 

(In  fact,  this  does  not  have  to  be  serious  per  se  as  the  example  3.2  B,  p.  73 
Ch.  IV  of  [1]  shows.)  The  real  difficulty  comes  because  the  travelling 
wave  equations  are  no  longer  in  a  form  that  makes  them  automatically  gradient- 
like  (as  they  must  always  be  in  the  case  of  one  locus  and  two  alleles  no 


-34- 


matter  what  the  form  of  the  gi  ^  is  ) .  ( Gradient-like  equations  are 

discussed  in  §6,  Ch.  I  of  [1]  ). 

Because  the  spatially  independent  equations  always  have  a  gradient  form, 
the  gradient-like  character  is  present  for  very  large  c  (cf.  Lemma  3,3  C) 
but  not  for  small  c  (cf.  Lemmas  3.3  A  and  3.3  B  where  the  gradient- like 
character  with  respect  to  H  made  the  treatment  very  simple  ) . 

Actually,  the  Laplacian  term  is  also  a  gradient  term  and  one  might 
expect  the  sum  of  two  gradients  to  be  a  gradient.  The  hitch  is  that  the 
metric  with  respect  to  which  the  spatially  independent  equations  cure  a 
gradient  system  is  not  the  usual  one  (but  rather  one  suited  to  the  simplex ) . 

In  [  14]  the  "non-gradient"  aspects  of  a  model  from  Ecology  are 
circumvented.  It  is  possible  that  similar  methods  might  be  used  in  special 
cases  of  the  equations  treated  here. 
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